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ABSTRACT 

The  objective  of  this  paper  is  to  study  the  basic  analytical  propagation  of  fundamental  and  higher  order  solitons  in 
the  optical  fiber.  The  characteristic  feature  of  soliton  formation  occurs  due  to  the  balance  between  Self-Phase  Modulation 
(SPM)  and  Group  Velocity  Dispersion  (GVD).  We  simulate  the  fundamental  and  higher  order  soliton  in  time  and 
frequency  domain  to  investigate  the  propagation  of  solitonic  pulses  in  nonlinear  optical  fibers. 
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1.  INTRODUCTION 

The  word  soliton  refers  to  special  kinds  of  wave  packets  that  can  propagate  undistorted  over  far  off  distances. 
Solitons  have  been  discovered  in  many  branches  of  Physics.  In  the  context  of  optical  fibers,  solitons  are  not  only  of 
fundamental  interest  but  also  have  practical  applications  in  the  field  of  fiber-optic  communications.  James  Scott  Russel 
was  the  first  scientist  who  observed  a  soliton  wave  in  1834  when  he  accidentally  noticed  in  the  narrow  water  canal  a 
smoothly  shaped  water  heap  that  to  his  surprise  was  able  to  propagate  in  the  canal  without  an  apparent  change  in  its  shape 
a  few  kilometers  along.  The  fact  of  propagation  of  this  solitary  wave  was  not  understood  for  a  long  time  until  appropriate 
mathematical  model  was  conceived  in  the  1960's  together  with  a  way  of  solving  nonlinear  equation  with  the  help  of 
inverse  scattering  method. 

The  phenomenon  of  modulation  instability  reveals  that  propagation  of  a  continuous-wave  (CW)  beam  inside 
optical  fibers  is  inherently  unstable  because  of  the  nonlinear  phenomenon  of  SPM  and  leads  to  formation  of  a  pulse  train  in 
the  anomalous  dispersion  regime  of  optical  fibers. 

2.  MATHEMATICAL  MODELING  AND  PROPAGATION  OF  THE  SOLITARY  WAVE  IN  OPTICAL 
FIBER 

The  propagation  of  light  can  be  described,  more  precisely,  with  Maxwell  mathematical  equations.  When  equations 
for  magnetic  and  electric  fields  are  combined  together  we  find: 

where  c  is  the  speed  of  light  in  the  vacuum  and  e0is  the  vacuum  permittivity.  The  induced  polarization  P  consists 
of  two  parts: 

P(r,t)=~PL(r,t)+P^L(f,t)  (2) 
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where  PL(r,t)  and  PNL(r,t)  are  related  to  electric  field  by  relations: 

K(r,  t)=  e0/_+J  j(l)(t-t')J(f,t')dt'  (3) 

iV(f,  0  =  £0/lC J(3)(t-t,,t-t2,t-t3)  *  E(?,U)E (f,h)E (?,t3)  dtidtadtj  (4) 
where  x(1)  and  /(3)  are  the  first  and  third  order  susceptibility  tensors. 

For  better  understanding  of  a  soliton  pulse  propagation  in  optical  fiber,  it  necessary  to  set  up  our  modelling  on  the 
mathematical  expression  (1).  We  suppose  that  a  solution  for  electric  filed  E  has  a  form  [1]: 

E(r,t)=A(Z,t)F(X,Y)exp(i(3oZ)  (5) 

where  F(X,Y)  is  transverse  field  distribution  that  corresponds  to  the  fundamental  mode  of  single  mode  fiber. 

A(Z,t)  is  along  propagation  axis  Z  and  on  time  t  dependent  amplitude  of  the  mode.  After  some  mathematical 
manipulations  we  can  come  to  the  equation  that  governs  pulse  propagation  in  optical  fibers  [1]: 

dA      _      dA     ip2  d2A      .   .  ,  ,2  . 

—  +Bi  — +  —  TT  =  iy\A  2A  (6) 

The  parameters  P[andp2  include  the  effect  of  dispersion  to  first  and  second  orders,  respectively.  Physically, 
P!=l/vg,  where  vg  is  group  velocity  associated  with  the  pulse  and  takes  into  account  the  dispersion  of  group  velocity. 
For  this  reason,  p2  is  called  the  group  velocity  dispersion  (GVD)  parameter.  Parameter  y  is  nonlinear  parameter  that  takes 
into  account  the  nonlinear  properties  of  a  fiber  medium.  Parameter  P^s,  in  real  case,  always  positive  but  on  the  other  hand 
parameters  XZD  and  y  can  be,  in  some  specific  case,  either  positive  or  negative.  The  parameter  P)is  closely  associated  in 
practice  with  better  known  parameter  called  dispersion  parameter  -  D  (ps/nm/km).  The  relation  between  them  is  in  the 
form  [1]: 

As  we  know,  dispersion  parameter  D  is  a  monotonically  increasing  function  of  wave-length,  crossing  a  zero  point 
at  wave-length  XZD,  which  is  called  a  zero  chromatic  dispersion  wave-length.  If  a  system  operates  with  wave-lengths 
above  XZD,  where  D  is  positive,  /?2must  be  negative  and  a  fiber  is  said  to  work  in  anomalous  dispersion  mode.  If  a  fiber  is 
operated  below  XZD,  the  D  is  negative  and  /?2must  be  positive.  In  this  case  a  fiber  is  said  to  operate  in  normal  dispersion 
mode.  As  regards  the  nonlinear  parameter  y,  it  can  generally  be  either  positive  or  negative,  depending  on  the  material  of  the 
wave  guide. 

For  silica  fiber  parameter  y  is  positive  but  for  some  other  materials  it  can  be  negative.  More  specifically, 
equation  (6)  has  only  two  solitons  -  either  dark  or  bright  solitons.  The  bright  soliton  corresponds  to  the  light  pulse  but  dark 
soliton  is  rather  a  pulse  shaped  dip  in  CW  light  background.  In  other  words,  the  dark  soliton  is  in  fact  negation  of  the 
bright  soliton.  While  there  is  maximum  of  light  in  the  bright  soliton,  the  dark  soliton  has  minimum  of  the  light. 
The  bright  soliton  can  propagate  in  only  such  a  waveguide  where  there  is  either  the  positive  nonlinearity  parameter  and 
fiber  in  anomalous  dispersion  regime  or  the  negative  nonlinear  parameter  but  in  normal  dispersion  regime. 
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3.  SIMULATION  AND  RESULTS 

Equation  (6)  can  be  defined  in  the  form: 

.du    s      d2u   .  2,  n 

l-   -T+  u  lu=0 

dz     2  St2 

Through  a  simple  conversion: 

t  =  (t-^ZVTo,  z=Z/LD,  u=V|y|Z,DA 


(8) 


(9) 


where  T0  is  pulse  width  and  LD=T0  /IP2I  is  the  dispersion  length.  Through  inverse  scattering  method  it  is  revealed 
that  the  solution  of  above  mentioned  equation  has  a  form: 


u(z,t)  =  N*  2/(eT  +  e_T)*  elz/2=  Nsech(T)eiz/2 

If  N  is  integer,  it  represents  the  order  of  the  soliton  pulse.  A  very  interesting  situation  comes  when  N=l. 


(10) 


(a)  (b) 
Figure  1:  Evolution  of  the  Pulse  Spectrum  for  N=l:  (a)  Time  Domain  and  (b)  Frequency  Domain 

In  the  case  of  first  order  soliton,  the  pulse  does  not  change  its  shape  at  all  as  it  propagates  in  optical  fiber. 
In  contrast  when  N  is  higher  than  one,  pulse  shape  is  not  stable  and  changes  periodically  with  soliton  period  Z0=  (n  /2)LD 
At  the  end  of  every  period  Z0  the  soliton  resembles  its  initial  simple  pulse  shape. 


(a)  (b) 
Figure  2:  Evolution  of  the  Pulse  Spectrum  for  N=3:  (a)  Time  Domain  and  (b)  Frequency  Domain 
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Figure  3:  Evolution  of  the  Pulse  Spectrum  for  N=4:  (a)  Time  Domain  and  (b)  Frequency  Domain 


It  is  evident  that  for  telecommunication  purposes  the  soliton  of  first  order  is  most  suitable  because  in  this 
application  it  is  necessary  to  keep  a  pulse  shape  stable.  N  defines  the  order  of  soliton  as  under: 


where  TO  corresponds  to  input  pulse  width,  P0  is  pulse  peak  power.  p2  takes  into  account  group  velocity 
dispersion  and  y  is  nonlinear  parameter  of  the  fiber  material. 


In  this  paper  a  numerical  analysis  based  on  simulations  has  been  employed  to  elaborate  the  propagation  of 
solitonic  pulses  in  nonlinear  optical  fibers.  The  Analysis  have  been  carried  out  in  time  and  frequency  domain  to  fully  study 
the  nonlinear  effect.  Several  discussions  based  on  numerical  results  have  been  conducted. 
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